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Solar electric propulsion (SEP) is being used for a variety of planetary missions sponsored by ESA, JAXA, and 
NASA and nuclear electric propulsion (NEP) is being considered for future, flagship-class interplanetary missions. 
Radioisotope electric propulsion (REP) has recently been shown to effectively complement SEP and NEP for 
missions to high-AU targets with modest payload requirements. This paper investigates the application of an 
advanced REP for a sample return from the comet Tempel 1. A set of mission and system parameters are varied 
with the goal of quantifying their impact on total mission payload. Mission parameters considered include trip-time 
and Earth return entry interface speed of the sample return system. System parameters considered include launch 
vehicle, power level of spacecraft at beginning of mission, and thruster specific impulse. For the baseline case of 
Atlas 401 and REP power level of 750 W, the mission time was 12 years, the payload was 144 kg, and the missions 
optimized to a single specific impulse generally within Hall ion thruster range. Other cases were investigated in 
support of graduate studies, and include the larger Atlas 551 launch vehicle and extended power level to 1 kW. The 
Atlas 551 cases tended to optimize dual specific impulses generally in the Hall ion thruster range for both legs of the 
mission. A power level of at least 1-kW and trip-time of approximately 11 years was required to obtain a total 
science payload close to 320 kg for the Atlas 401 launch vehicle. An Atlas 551 launch vehicle yielded a science 
payload of approximately 540 kg for the case of 1-kW of power and an 11 -year trip time, and nearly 250 kg of 
science payload for the case of 1-kW of power and a 6-year trip time. Results are also reported indicating the 
performance ramifications of meeting a reduced Earth entry interface velocity constraint. 


I. Introduction 


E lectric propulsion (EP) is being used increasingly for planetary missions. Solar electric propulsion (SEP) has 
recently been successfully demonstrated by NASA’s Deep Space l 1 , ESA’s SMART l 2 , and JAXA’s 
HAYABUSA 3 spacecraft. Nuclear electric propulsion (NEP) is also being considered for missions requiring large 
payload masses and powers. SEP and NEP have been shown to provide major values for a range of planetary 
missions. In addition, radioisotope electric propulsion (REP) has been shown to be a strong candidate for planetary 
missions with small to modest payloads traveling to high AU. Early studies by Noble 4,5,6 and subsequent NASA 7,8,9 
work analyzed REP efficacy to the outer solar system and near-interstellar space. Advantages of REP include a 
power source that is basically insensitive to AU, enables the accomplishment of modest missions in reasonable 
mission times, and has a legacy of successful applications on planetary missions. 

This paper investigates the feasibility of using REP for a sample return mission to comet Tempel 1 in the 2010 
time frame. For this analysis, two-body dynamics govern the orbit transfer, and the propulsion system is optimized 
for parameters such as launch energy ( C 3 ) and the REP specific impulse (/sp). Several mission scenarios that 
maximize the mass returned to Earth for a range of initial REP power levels and sample return earth entry speeds are 
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obtained. Using systems analysis, the spacecraft size, propulsion subsystems, and net payload are obtained for the 
various mission cases considered. The baseline case shown in this paper represents early results in support of a 
study performed for NASA’s Science Mission Directorate’s Radioisotope Electric Propulsion Program 10 . This 
paper also reports later results obtained through support of graduate work at the University of Missouri. 


II. System and Mission Analyses 

A. System Approaches 

The overall spacecraft model is based on state-of-art vehicles such as Mercury MESSENGER and DAWN. The 
REP power system (RPS) consists of a set of radioisotope-based general-purpose heat sources, a Stirling energy 
conversion system that transforms the heat to electrical power, power conditioning components, and distribution 
cabling. The RPS power source was assumed to be a second- generation concept with a specific power of 8 W/kg 
that provided power to the electric propulsion system and a fixed 100 W to the spacecraft. For all calculations, it is 
assumed that the REP power degrades at a constant rate of 1 . 15% per year. 

The electric propulsion system consists of a power-processing unit that converts the bus power to the required 
thruster operating currents and voltages, the electric thrusters, and cable harnesses from the power processors to the 
electric thrusters. The electric propulsion hardware masses were provided by the Glenn Research Center. The system 
includes two PPUs, one being a spare, with fixed specific masses of 12 kg/kW. The thruster specific masses are 3.25 
kg/kW, 3.00 kg/kW, and 2.75 kg/kW for 500 W, 750 W, and 1000 W respectively. The thruster is an advanced 
design assumed to have an operating life of 30,000 hrs. The thruster specific impulse is limited to a range of 1000 to 
4000 s, and a simple parametric relation is used to compute thruster efficiency (rf) as a function of 7 sp : 
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This parametric relationship was developed by the Glenn Research Center specifically for this study. Thrust 
magnitude (7) is computed as a function of input power (P), efficiency ( rj ), and specific impulse: 


r = 


2 tjP 
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This range of covers both the theoretical ranges of a Hall thruster (defined for this study as 1000 s to 2500 s) 
and of a gridded ion thruster (defined for this study as 2500 s to 4000 s). An actual thruster cannot operate over the 
entire / sp range selected. Trajectories that have an outbound 7 sp in one range and an inbound 7 sp in another range do 
not include a system mass penalty, which would be needed to reflect the realities of having two separate propulsion 
systems. It is assume that the REP power degrades at a constant rate of 1.15% loss per year. 

Injected mass for the Atlas 401 and 551 launch vehicles is computed as a function of launch energy, as shown 
in Fig. 1. 




B, Mission Profile and Assumptions 

The interplanetary trajectory begins with injection into heliocentric space by either an Atlas 401 (baseline 
mission) or Atlas 551 launch vehicle, providing the spacecraft with the desired launch energy (C 3 ). Figure 1 shows 
the capabilities of both launch vehicles. Immediately after leaving geocentric space, the ion propulsion system is 
started and a powered arc is initiated. It is assumed that the spacecraft is equipped with multiple xenon ion engines 
that have performance and lifetime characteristics similar to a small advanced Hall or gridded ion thruster. Because 
the REP system can only supply 500-1000 W of power for the EP system, it is assumed that only a single ion 
thruster is operated at a time. After a prolonged bum, erosion eventually renders the working thruster unusable, and 
the next thruster is started to continue the powered arc. Therefore, thruster lifetime issues or maximum propellant 
limitations do not impose additional constraints on the trajectory design as long as total engine mass is properly 
accounted for. 

After a continuously powered outbound cruise phase, the spacecraft matches the heliocentric orbit of comet 
Tempel 1 and spirals down to a stable orbit about the comet. A small probe (maneuvering sampling system) 
descends to the comet surface and returns samples to the orbiting primary spacecraft. It is assumed that the main 
spacecraft leaves 50 kg of the sampling system at the comet after sampling operations. A 60-day stay-time 
constraint at Tempel 1 is imposed in order to account for rendezvous, science, and sampling operations. After a 
minimum stay time of 60 days, the Earth-return leg is initiated and thrusting is maintained until the spacecraft 
reaches proper orbital conditions for a ballistic coasting trajectory to Earth. The sample return mission ends with the 
sample return system’s direct ballistic entry into Earth’s atmosphere. The primary spacecraft is assumed to be 
diverted from an intercept with Earth into a heliocentric orbit. 

C. Trajectory Optimization 

The comet sample return mission design objective is to maximize the science payload mass. The trajectory 
design variables include launch date, injection point into optimal Earth departure condition, outbound and return 
powered arc durations, thrust direction during the powered arcs, specific impulse for the outbound and return legs, 
Tempel 1 rendezvous date, stay time at the comet, and Earth-return date. The powered equations of motion are 




numerically integrated by using a fixed-step, fourth-order Runge-Kutta routine in a heliocentric frame. Coast arcs 
are computed analytically. Therefore, the outbound and inbound trajectories begin and end at the respective body’s 
sphere of influence. Details of the numerical simulation of the interplanetary trajectory are presented in Ref. 10. 

Numerical solutions of the maximum-delivered-mass problem are obtained by using a direct trajectory 
optimization method (DTOM) 11 . The optimal control problem is replaced by a nonlinear programming problem, 
which is solved by using sequential quadratic programming (SQP), a constrained parameter optimization method. 
The SQP code used here is described in Ref. 12 and uses forward finite differences to compute gradient information. 
The thrust direction is defined by three direction cosines relative to a rotating radial-transverse-normal frame and the 
time histories of these direction cosines are parameterized by linear interpolation through a set of discrete nodes. 
Midpoint rendezvous conditions are imposed by six SQP equality constraints that maintain a match between the 
orbital elements of the spacecraft and Tempel 1 at the rendezvous date. Three additional SQP equality constraints 
maintain a heliocentric position match between the spacecraft and Earth at the terminal end of the trajectory 
(ballistic entry). Finally, a single inequality constraint is imposed on the hyperbolic approach speed at the Earth’s 
sphere of influence in order to maintain feasible limits on the entry speed. The orbital elements of the Earth and 
Tempel 1 are computed using an accurate solar system ephemeris. 13 

Ill, Analyses Results 

A. Payload Variation with Varying Initial Power 

Our initial trials involve computing optimal sample return missions over a range of initial power levels for the 
REP system. The initial beginning-of-life power is set to 500, 750, and 1000 W and the optimal trajectories that 
maximize science payload mass at Earth return are obtained. For these initial cases, it is assumed that the hyperbolic 
approach speed for the ballistic Earth-return leg is constrained to <10.7 km/s, which corresponds to 15.4 km/s at 
entry interface (125 km above the Earth’s surface). In a later section of this report, results are shown that indicate 
the mission impact of reducing the Earth entry interface speed to <13 km/s. In addition, the launch energy is fixed 
to obtain optimal missions over a fixed set of C 3 . 

Figure 2 shows a representative mission profile for the case using the Atlas 401, where initial power is 750 W 
and C 3 = 36 km 2 /s 2 . Because the REP power is so low, the spacecraft requires over two revolutions about the sun 
(with continuous thrust) before it sufficiently increases energy to reach Tempel L After the 60-day stay time at 
Tempel 1, a single powered arc (one-half revolution) decreases energy and the spacecraft begins the transfer back to 
Earth. The final phase of the mission is a ballistic coasting trajectory that ends with an Earth atmospheric entry 
speed of 15.4 km/s. The total trip time is about 1 1 years (as compared to 12 years for the baseline mission), and this 
mission profile and relative geometry appear to provide a good local minimum for the launch year of 2010. 

Figure 3 presents the maximum spacecraft science payload as a function of launch energy for the three initial 
power levels of the Atlas 401. It is clearly seen that an increase in power level for the REP spacecraft results in 
greater science payload mass. For each power level, there exists an optimal launch energy, and the optimal C 3 tends 
to decrease as the power increases. These trends in optimal final mass and optimal launch energy are easily 
explained: higher power increases thrust magnitude, and therefore the REP system can perform a greater share of the 
required orbit transfer, which subsequently reduces the magnitude of the energy contribution from the launch 
vehicle. Therefore, the decrease in launch energy increases the injected mass, which ultimately increases the 
payload at Earth return. 



Figure 2. Representative trajectory for Tempel 1 sample return mission: initial power of 750 W, Atlas 401 

with C 3 =3 6 km 2 /s 2 



Figure 3. Maximum payload mass for different initial power levels and launch energies (Atlas 401), 




Table 2 summarizes the best mission designs for the Atlas 401 for each initial power level and provides the 
mission characteristics for the three values of C 3 that give the best final mass solution for each power level shown in 
Fig. 5. A distinct trend among the different missions became apparent: all missions (regardless of initial power) use 
very nearly the same mission profile with the same launch date, rendezvous date. Earth-arrival date, and 11-yr trip 
time. Therefore, relative geometry between Earth and Tempel 1 in the 2010 opportunity dictates the best trajectory. 
Furthermore, when power level is increased, the REP system is able to operate at a higher specific impulse, which 
increases the overall efficiency of the transfer and reduces the xenon propellant requirement. For example, note that 
when initial power is 500 W, the mission requires an / sp of 1267 s, which is close to the lower bound of 1000 s. This 
low specific impulse is required to increase thrust in order to reach Tempel 1 with a continuous powered arc and the 
optimal flight profile shown in Fig. 2 (i.e., two and one-quarter revolutions about the sun for the optimal rendezvous 
position with Tempel 1). As power is increased, optimal J sp settings increase. 


Table 2. Mission parameters for optimal sample return missions using the Atlas 401 

Earth Entry Speed = 15.4 km/s 


Initial 

Power 

(W) 

C 3 

(kmVs 2 ) 

Launch 

Date 

Trip Time 
(yrs) 

mg 

(kg) 

®prop 

(kg) 

Payload 

(kg) 

^sp 

(s) 

500 

42.6 

5/28/2010 

11.01 

1359 

614 

158 

1267 

750 

36.9 

5/23/2010 

11.03 

1572 

577 

259 

1773 

1000 

34.3 

5/20/2010 

11.04 

1675 

514 

316 

2316 


Next, optimal missions were determined using the Atlas 551 launch vehicle. Again, hyperbolic approach speed 
at Earth was constrained to 10.7 km/s (15.4 km/s at Earth entry interface). Table 3 summarizes these best mission 
designs for the optimal launch energy for a range of initial power levels. The larger Atlas 551 vehicle permits 
different options for rendezvous dates and trip times. For example, the larger launch vehicle allows for cases with 
significantly higher launch energy than the Atlas 401. The higher launch energy effectively eliminates one 
revolution about the sun on the outbound leg, and subsequently the trip time can be reduced by 3 or 6 years for the 
each of the cases. Ail missions have an Earth-return date of early June in the years 2022, 2019, and 2016. Figure 4 
shows a trajectory example for the previously mentioned 6 year direct outbound transfer to Tempel 1 for the case 
where initial power is 1000 W. Note that for the outbound and return legs, the total transfer angle is around 360 deg. 
The science payload mass for this high-power case is only modestly less than the corresponding Atlas 401 case, but 
the trip time has been cut nearly in half from 11 years to 6 years. Table 3 also shows an interesting case (power 
level is 1000 W) where the trip time is 11 years, and the payload has reached 535 kg, which is almost double the 
payload of the corresponding Atlas 401 case. 

Table 3. Mission parameters for optimal sample return missions using the Atlas 551 

Earth Entry Speed = 15.4 km/s 


Initial 

Power 

(W) 

C 3 

(km 2 /s 2 ) 

Launch 

Date 

Trip Time 
(yrs) 

mo 

(kg) 

Ulprop 

(kg) 

Payload 

(kg) 

'sp 

(outbound/ 
return legs) 
(s) 

500 

70 

5/31/2010 

11.01 

1612 

780 

184 

1000/1834 

750 

48 

6/3/2010 

11.00 

2575 

1261 

361 

1000/1967 

1000 

40 

5/31/2010 

11.01 

3057 

1330 

535 

1243/2422 

500 

79 

5/10/2010 

9.07 

1260 

692 

51 

1000/1089 

750 

66 

5/2/2010 

9.09 

1745 

887 

151 

1232/1212 

1000 

64 

4/30/2010 

9.09 

1850 

764 

257 

1739/1592 

500 

85 

5/9/2010 

6.08 

1054 

526 

54 

1000/1189 

750 

68 

5/1/2010 

6.10 

1678 

858 

142 

1003/1234 

1000 

61 

4/24/2010 

6.12 

1953 

869 

247 

1300/1491 



Payload, kg 



X, AU 

Figure 4. Six-year Tempel 1 sample return mission: initial power of 1000 W, Atlas 551 with C 3 =61 km 2 /s 2 
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Figure 5. Trip Time versus payload trade for Atlas 551 and Atlas 401 while varying power levels 




The data of Tables 1 and 2 are now organized into a payload versus trip time plot shown in Figure 5. A 
comparison of the Atlas 401 and 551 launch vehicles shows that for the 55 1 launch vehicle, a significant decrease in 
trip time can be realized over the Atlas 401 cases. Likewise, for the same trip time, the 551 can deliver more science 
payload mass than the 401. Yet in the case of 500 W REP, the difference in science payload delivery capability is 
minimal. For the 1 1 year trip time, the 750 W REP spacecraft and the 551 launch vehicle can provide significantly 
more payload than the corresponding case for the 401 launch vehicle. The 1000 W spacecraft injected by the 551 
launch vehicle can deliver about the same payload as the 1000 W Atlas 401 case, but at a very significant 6 year 
reduction in trip time. If large payloads (on the order of 500 kg) are required, then a longer 1 1 year trip time and the 
maximum 1000 W of power will be needed. 

B. Payload Variations for Earth Return Entry Interface Speed 

Next, several optimal sample return missions were obtained for a range of constraints on hyperbolic approach 
speed at Earth return. Table 4 shows a parametric study on reducing entry speed at Earth interface for the Atlas 401 
with initial REP power of 1000 W. These 11-yr mission cases will allow us to observe the trends in payload mass 
when Earth entry speed is varied. Note that entry speed can be reduced by adding trip time for another revolution 
about the sun, resulting in an energy reduction. Figure 6 shows a representative 13-yr mission that reduces the 
reentry speed to 12.5 km/s, which is closer to the reentry speed of the Apollo capsule. The baseline entry speed of 
15.4 km/s is aggressive and may not be achievable by 2010, Note that an additional revolution around the sun is 
needed for this additional energy depletion before reentry. The payload for this 13-year case is, however, greatly 
reduced to an inadequate 53 kg. Table 4 shows that the final mass dramatically decreases as the limit on reentry 
speed is reduced. 

Table 4. Optimal sample return missions using the Atlas 401, initial power level of 1000 W, and reduced 


Entry 
Speed, V E i 
(km/s) 

c 3 

(km 2 /s 2 ) 

Launch 

Date 

Trip 

Time 

(yrs) 

mo 

(kg) 

Ulprop 

(kg) 

Payload 

(kg) 

15.0 

34 

5/19/2010 

11.04 

1673 

514 

314 

14.5 

35 

5/19/2010 

11.04 

1664 

518 

307 

14.0 

35 

5/20/2010 

11.02 

1657 

535 

291 

13.5 

35 

5/22/2010 

11.04 

1636 

583 

247 

13.0 

35 

6/7/2010 

12.96 

1668 

826 

100 

12.5 

35 

6/8/2010 

12.97 

1623 

846 

53 


Y, AU 


4 
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Figure 7. 13-yr Tempel 1 sample return mission with 12.5 km/s Earth entry interface speed and Earth 

gravity assist during return leg. 




Figure 7 shows a trajectory option that uses an Earth gravity assist (EGA) during the return leg in order to 
further reduce the Earth entry interface speed. An EGA may or may not be feasible with the use of a nuclear power 
system, but this option is being presented for consideration and comparison. Note that the EGA follows a ballistic 
return arc, and that a relatively short powered arc is required after the EGA for proper entry conditions. It may be 
possible to jettison the REP system before the first EGA, and then use multiple phased EGAs to eventually set up 
the desired entry conditions. 

Table 5 shows the mass benefits of including one EGA before reentry. Note that payload mass is improved 
when compared to the Atlas 401 cases without a gravity assist (Table 4), with the expense of additional trip time (the 
mission is now 13 years), and that solutions are possible with 750 and 1000 W of power. 


Table 5, Optimal sample return missions using the Atlas 401 and Earth gravity assist for further reduced 

Earth entry interface speed 


Initial 

Power 

(W) 

Entry 
Speed, Vn 
(km/s) 

c 3 

(km 2 /s 2 ) 

Launch 

Date 

Trip Time 

(yrs) 

m« 

(kg) 

m prop 

(kg) 

Payload 

(kg) 

750 

12.5 

38 

5/23/2010 

12.87 

1549 

613 

217 

1000 

12.5 

35 

5/20/2010 

12.97 

1654 

547 

277 

1000 

13.0 

35 

5/20/2010 

13.00 

1655 

541 

282 


Lastly, Table 6 presents 6-yr mission solutions using the Atlas 551 for a range of reduced Earth entry interface 
speed. Only for the 1000 W power case was a reduced entry speed solution found. It is only possible to reduce 
entry speed to 14 km/s using the 6-yr mission profile with the Atlas 551. 

Table 6. Optimal sample return missions using the Atlas 551, initial power level of 1000 W, 6-yr mission 

time, and reduced Vn 


Entry 
Speed, Vex 
( km/s) 

c 3 

(km 2 /s 2 ) 

Launch 

Date 

Trip 

Time 

(yrs) 

mo 

(kg) 

Mprop 

(kg) 

Payload 

(kg) 

15.5 

61 

4/24/2010 

6.12 

1953 

869 

247 

15.0 

61 

4/24/2010 

6.13 

1959 

907 

225 

14.5 

61 

5/20/2010 

6.13 

1972 

1000 

170 

14.0 

61 

5/22/2010 

6.14 

1967 

1114 

91 


TV. Conclusions 

Several designs for a comet sample return mission using an advanced REP based ion propulsion system have 
been obtained. A set of mission and system parameters were varied with the goal of quantifying their impact on 
total mission payload. Mission parameters considered include trip-time and Earth return entry interface speed of the 
sample return system. System parameters varied include launch vehicle, power level of spacecraft at beginning of 
mission, and thruster specific impulse. 

For the baseline 750 W case, the specific impulse of the EP thruster was optimized so that the trajectories 
exhibit an approximately 12 year trip time for the Atlas 401 launch vehicle. An initial quick-look did not show 
significant improvements for the higher power level of 1-kW. Work in addition to the baseline work was later 
performed in support of graduate studies at the University of Missouri. In this additional work, improvement was 
found in transfer time to approximately 1 1 -yr total flight time for the same Atlas 401 launch vehicle for each of the 3 
power levels investigated. Further, cases were found providing a maximum science payload mass for a launch 
energy of C 3 = 33 km 2 /s 2 , C 3 = 36 km 2 /s 2 , C 3 = 41 km 2 /s 2 for system power level of 1000 W, 750 W, and 500 W 
respectively. A power level of at least 1000 W and trip- time of approximate 11 years was required to obtain a total 
science payload of approximately 280 kg for the smaller Atlas 401 launch vehicle. The larger Atlas 551 launch 
vehicle yielded a science payload increase to approximately 530 kg for the case of 1000 W of power and a 12 year 
trip time, and nearly 250 kg of total payload for the case of 1000 W of power and a 6 year trip time. For the Atlas 
401 launch vehicle cases, the outbound and inbound legs of the missions optimized to a single specific impulse 


generally in the Hall ion thruster range. The Atlas 551 cases tended to optimize generally in the Hall ion thruster 
range for both outbound and inbound legs of the mission. 
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